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ABSTRACT 

If the orientations of galaxies are correlated with large-scale structure, then anisotropic se- 
lection effects such as preferential selection of face-on disc galaxies can contaminate large 
scale structure observables. Here we consider the effect on the galaxy bispectrum, which has 
attracted interest as a way to break the degeneracy between galaxy bias and the amplitude 
of matter fluctuations erg. We consider two models of intrinsic galaxy alignments: one where 
the probability distribution for the galaxy's orientation contains a term linear in the local tidal 
field, appropriate for elliptical galaxies; and one with a term quadratic in the local tidal field, 
which may be applicable to disc galaxies. We compute the correction to the redshift-space 
bispectrum in the quasilinear regime, and then focus on its effects on parameter constraints 
from the transverse bispectrum, i.e. using triangles (^1,^2,^3) in the plane of the sky. We 
show that in the linear alignment model, intrinsic alignments result in an error in the galaxy 
bias parameters, but do not affect the inferred value of erg. In contrast, the quadratic alignment 
model results in a systematic error in both the bias parameters and erg. However, the quadratic 
alignment effect has a unique configuration dependence that should enable it to be removed 
in upcoming surveys. 

Key words: large-scale structure of Universe - cosmology: theory. 



1 INTRODUCTION 

While the evolution of dark matter perturbations in the current 
ACDM model is well understood theoretically, the relation between 
the galaxy distribution and the large scale (dark) matter distribution 
is complicated by the detailed physics of galaxy formation and dif- 
ferent models may lead to different clustering properties of galax- 
ies. In particular, while local theories of galaxy formation predict 
that the galaxy density fluctuations trace the matter fluctuations on 
large scales, they also predict that the two are related by the bias 
parameter b, which is in general not known a priori ( Kaiser 19840. 
The unknown bias parameter represents a key problem for attempts 
to measure the growth of cosmological perturbations using galax- 
ies. 

In combination with the galaxy power spectrum, third or- 
der galaxy clustering measures such as the bispectrum or (equiv- 
alently) the 3-point correlation function can be used to mea- 
sure non-linear galaxy bias and break the degeneracy between 
the normalization of the matter power spectrum, trg, and the 
linear galaxy bias. This enables one remove the effects of 
galaxy biasing and measure the c osmological growth of struc- 
ture from the galaxy distribution (jFrvl Il994l : IVerde et al] 1 19981 : 
IScoccimarro et al. 1999; Verde et al. 2000), and thus constrain dark 
energy (e.g. lDoln ev et al. 2006). Recently third order galaxy clus- 
tering has_been_analj^d_by sever^ the bispec- 
trum iScoccim arro et al.l I2OOII : iFeldman et"ZI I2OO1I : IVerde et al] 



jjing & B6me3l2004 iKavo et aljliooj : iNichol et alj|2006l) Usin 



l2002l : lKulkariii et al. 200'^ and the three point correlation function 



mock catalogs from numerical simulations, ISefusatti et alj j200i 
show that a combined analysis of the galaxy power spectrum and 
bispectrum including their cross-correlation contains significant in- 
formation on galaxy bias and fundamental cosmological parame- 
ters and helps break parameter degeneracies of other cosmological 
probes. 

The most important systematic errors in interpreting the ob- 
served galaxy clustering arise in the non-linear regime,where the 
behavior of galaxy biasing and models of the (redshift space) 
galaxy power spe ctrum and bispectrum are difficult to model (see 
ISmith et^l2008l for the complicati ons of a curre nt model of the 
redshift space Bispectrum). Recentlv fairatal j2009h showed that the 
alignment of galaxies by large scale tidal fields can cause a system- 
atic error in the determination of the linear redshift space distor- 
tion parameter yS tKaiser..l987) : the alignment of galaxies with the 
tidal field (along the stretching axis of the field for large ellipti- 
cal galaxies) in combination with a viewing direction dependent 
galaxy selection effect, e.g. preferential selection of galaxies which 
are observed along their long axis, will lead to a selection probabil- 
ity for galaxies which is modulated by the tidal field along the line 
of sight. This results in an anisotropy in redshift-space clustering 
with the same scale and angular dependence as the linear redshift- 
space effect. In this paper we will explore the implications of such a 
tidal alignment contamination for the observed galaxy bispectrum 
and how it affects the measurement of galaxy bias parameters. 

Throughout this paper we assume a standard ACDM cosmol- 
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ogy with the best-fit WMAP 7 jKomatsu et aljlloioh parameters, 
and assume Gaussian initial density perturbations. 

We begin in Section |2] with a derivation of the standard red- 
shift space galaxy bispectrum and discuss toy models for physical 
processes that cause alignments of galaxy orientations with large 
scale structure. In Section [3] we explain how tidal alignments of 
galaxies in combination with an orientation dependent galaxy se- 
lection modify the observed galaxy distribution and calculate the 
corresponding corrections to the galaxy bispectrum. Using a Fisher 
matrix technique we then estimate the systematic error induced by 
tidal alignments to measurements of galaxy bias parameters from 
angular clustering in Section^ We conclude and discuss mitigation 
strategies in Section|5] 



2 THEORETICAL BACKGROUND 

In this section we derive the redshift space galaxy bispectrum 
to second order in pe rturbation theory (for a review, see e.g. 
ISemardeau et al.ll2002[) . and discuss toy models for the aligimient 
of galaxies with the large scale tidal field. 



2.1 Galaxy bispectrum 

The matter bispectrum B is defined as 

{6{ki)6{k2)~6(ki)) = {iTtfSo (A:,23) 5(A:,, A:,, k^) , (1) 

where S{k) is the matter density contrast in Fourier space, 6\) the 
Dirac delta function, and ki^i = ki + k2 + kj. The bispectrum 
vanishes for a Gaussian random field. 

To second order perturbation theory the density contrast is 
given by 



S(k) = S^^\k) + 



J {2ny 



F2(k,,k-kt)~S'-'\kiW»(k-k,), (2) 



with 6^^\k) the linear density contrast, and the second order density 
kernel 



Fiikuki) = - + —— 
1 2k[k2 



h kA 2 I kj k: 
k2^ h]^ l\ hk2 



(3) 



Hence the matter bispectrum induced by non-linear gravitational 
evolution at tree-level is given by 



B{k^,k2,ki) = 2F2(ki, k2)P(ki)P(k2) + 2 perm.. 



(4) 



where P(k) is the linear matter power spectrum, k^ = —k\ - and 
"2 perm." indicates that the 2 permutations (^2, ^3) and (Ai, A3) are 

also included in the summation. 

Using the local bias approximation (e.g. IFfv & Gaztafiagal 

Il993h . the galaxy density contrast Sg can be expressed as a non- 
linear function of the matter density contrast 



5„{x) = bi6(x) + ]^b26{xf + ■ 



(5) 



where the expansion coefficients are the linear {bi) and non-linear 
galaxy bias factors. In reality, galaxy biasing may be mor e com- 
plicat ed, especially on small scales, due to I-halo terms l lSeliald 
l200d ) and nonlocal dependences such as the stren gth of the lo- 
cal tidal field ( McDonald 200^: iMcDonald & Rovii2009 ). How- 
ever, in simulations the local bias model is found to be a fair de- 
scription of non-linear halo clustering on large scales with an ac- 
curacy of a few percent (e.g. lMarfn et al.ll2008l : IGuo & Jing||2009l : 



anera & GaztanagalboogI) . which is sufficient at the level of this 
analysis. 

Then the galaxy bispectrum is related to the matter bispec- 
trum via 

Bg{ki,k2, ki) ^ b\B{ki,k2, k^) + b\b2 {P{ki)P{k2) + 2 perm.) , 

(6) 

and similarly for the galaxy power spectrum, 

Pg(k) = b\P{k). (7) 

To arrive at an expression for the redshift space galaxy bis- 
pectrum we have to transform radial coordinates to redshifts space. 
In the plane-parallel approximation, the mapping from real space 
position X to coordinate jc" in redshift space is given by 



h ■ u(x) , 

X H fl. 

Ha 



(8) 



where u(x) is the peculiar velocity field, and h is the direction of 
the line of sight. The velocity field is curl-free, V x u(x) = 0, at all 
orders in perturbation theory. Its divergence is given to linear order 
in perturbation theory by 

ik ■ u"\k) = aHf6^'\k), (9) 

where / = dln(G)/dln(a) is the logarithmic growth rate of linear 
perturbations (equal to roug hly Ci!^^ in general rela tivity). Higher- 
order contributions to V • « jBemardeau et al.ll2002h are analogous 
to Eq. lO, e.g. 

ik-u^^\k) = aHf 

-'J (27r)3 

with the kernel 



G2(ki,k-ki)d'^'>(ki)6^'>{k-kt), (10) 



ko 



G2{k,,k2)=l + ^— 
1 2k\k2 



k\ fe 



4 
-I- - 

7 



k\- k2\ 



(11) 



Taking into account the Jacobian of this mapping of jc — > jc'* 
(Eq. (8)1, and approximating the peculiar velocity field by the sec- 
ond order bulk veloc i ty field, the galaxy densi ty is redshift space is 
( iHeavens et alJl998l : fScoccimarro et al]ll999l) 

%{V) = {bi+ffiWm 
d'k'i 



f ^Z2 (k],k' - k\yd^'Kk])d^\k' - k\). 



(12) 



where k" denotes a Fourier mode in red shift space, and fi = k ■ h 
is the cosine of the angle between the wave vector and the line of 
sight (we may analogously define fii , fii2, etc.). The mode-coupling 
function Zt is 



Z2(A:i, ^2) 



biF2(ki,k2) + fn]2G2(kuk2) 



tp.ih^+f^il)+f^(b,+ffi]) 
ki k2 



b2 

+ f (13) 



Hence we can write the redshift space galaxy bispectrum as 

Bl(k\,kl,kl) = 2(hi+f^])(h,+f^l)P(kl)P(kl)Z2(k\,kl) 

+2 perm. (14) 

Note that this expression does not include the Finger of God effect 
due t o the virialized motion of galaxies within a cluster (Jackson] 
Il972h . which is important when one of the kj has a large line-of- 
sight component. While this effect is important even on weakly 
non-lin ear scales, it is usually handled by phenomenological mod- 
els (e.g. Hatton & Cole'l998l: IVerde et al.lll998l : IScoccimarro et al.l 
J999 ; Peacock^et al. 2001), a compressio n of radial coordinates for 
galaxies living in the same cluster (e.g. iTegmark et aL|l2004h, or 
by re constructing the redshift-space halo density field ( Reid et al] 
l2009b . 
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2.2 Toy models of tidal alignments 

2.2.1 Halo shape distortion.^: linear alignment 

In the linear alignment model jCatelan et al.ll200ll) the shape and 
orientation of a galaxy are assumed to be determined by the shape 
of the halo it resides in. It is thought that the gravitational collapse 
of an initially spherical overdensity in a constant gravitational field 
leads to triaxial haloes, such that the halo will be prolate if the over- 
density is stretched by the large scale tidal field and oblate if it is 
compressed. This mechanism is believed to lead to a net correla- 
tion of halo orientations even though overdensities typically are 
not spherical, and such an alignment has been confirmed by sim- 
ulations (e.g. Faltenbacher et al. 2009). 

The relation between halo shape and galaxy shape is com- 
plic ated by galaxy formatio n and differs between galaxy types 
(e.g. lFaltenbacher et alj|2007l) . but at least for luminous red galax- 
ies (LRGs) there is observational evidence fo r an aligrmient 
of the LRG with the ma|or axis of its host teinggelj [1983 : 



iFaltenbacher et al.ll2007l ; iNiederste-Ostholt et alj 1 20 1 on. The re are 
also correlations with large-scale structure teinggelilfigsi ): with 
the Sloan Digital Sky Survey (SDSS) it has even been possible to 
measure the scale dependence of these correlations and show the 
consistency of their spectral index with the predictions of the linear 
tidal a ligrmient model and the ACDM power spectrum jHirata et al.l 
l2007h . 



2.2.2 Tidal torques: quadratic alignment 

The orientation of a disc galaxy is determined by the direction 
of its angular momentum, which builds up due to tidal torquing 
during early stages of galaxy formation if the proto-galaxy's iner- 
tia tensor is anisotropic and misaligned with the local shear field 
(lHovl^ll949U Sciami^'l955': 'Peeblej'1969': 'Poroshkevichl Il97d : 
I984lc'rit 



ite 



ICrittenden et al. 2001). See Schafer (2009) for a re- 
view of tidal torquing and the build up of angular momentum cor- 
relations. 

Following ITee & Perj j2000l) . we parameterize the correlation 
between moment of inertia and the shear field by 



\ + a 



-6ij - aTii,Ti,j 



(15) 



which is also the most general quadratic form possible. Here tjj 
is the unit normalised traceless tidal field tensor (TjjTjj = I) and 
a is a dimensionless coupling parameter, e.g. a = | at leading 
order in perturbation theory if shear and inertia tensor are mutually 
uncorrelated. It is also possible for a to be much smaller, e.g. if the 
angular momentum vecto r of the disk is on ly partially aligned with 
that of the host halo (e.g. van den Bosch et all 2 002). 

Note that in non-linear theory spin-induced alignments also 
have a linear contribution at large scales because the large- 
scale tidal field induces correlations of the small-scale tidal field 
and inertia tensor that lead to a nonzero contribution to (LjLj) 
jHui & Zhang|l20()8l) . although this linear effect has n ot been ob- 
served for late-type galaxies despite several sear ches jLee & Peril 
l2007l : lHirata et alj2007l : lMandelbaum et alj2009l) . 



3 TIDAL ALIGNMENT CONTAMINATION 

As discussed in the previous section, the orientation of galaxies 
likely is not random but correlated with large scale structure, and 
in combination with observational galaxy selection criteria which 



depend on the galaxy orientation relative to the line of sig ht, this 
may m odify the observable galaxy distribution. Following iHiratal 
(|2009|), we will now introduce the basic notation needed to dis- 
cuss galaxy orientation and viewing direction dependent selection 
effects. 

Let the galaxy orientation be described by the Euler angles 
(6, <p, if/) through a rotation matrix Q{d, <p, tf/). This matrix trans- 
forms "lab" frame coordinates to a coordinate system aligned 
with the galaxy. Due to tidal alignments the probability distribu- 
tion p{Q\x) for the orientation of a galaxy at position x may be 
anisotropic and a function of the local environment of x. The ob- 
servational galaxy selection probability depends on the direction of 
the line of sight, h, and the galaxy orientation, specifically on the 
direction of the line of sight in the galaxy frame Qh. We define 



Poc 1 + r(Qh,x), 



(16) 



where the anisotropic part T is zero when averaged over all possible 
galaxy orientations or viewing directions. 

The observable galaxy distribution A'( selected) hence is mod- 
ified compared to the true galaxy distribution A'(true) by 

A'( selected) 

-('"■ 



A'(true) 



h\x) oc r p(QW[l + T(Qn,x)] d^'Q 

Js0(3) 

= 1 + r p{Q\x)r(Qh,x) d^Q 

Js0(3) 

= l + e(nW, (17) 

which is the average of Eq. l ll6t over the distribution of galaxy 
orientations, and where we have defined the orientation dependent 
selection function E(h\x) in the last step. As the average of T over 
all galaxy orientations vanishes, Eq. dlVI ) implies that 6 vanishes 
if either the galaxy orientations are isotropically distributed or if 
the probability for selecting a galaxy is independent of Qh, i.e. if 
T = 0. 

The observed galaxy density is modified by the orientation 
dependent selection function such that 

1 + Sf\x'') = ([l + 6,(x)] [I + €mx)]J , (18) 

where the term in round brackets is the orientation modulated real 
space density of selected galaxies, and where the superscript s de- 
notes the transform to redshift space. Expanding to second order in 
the matter density field, this implies: 



5f >(F) + ?'<"(«!*') + Sf^\k') + e''^\h\k') 



<5J*"(A5)r*"(«|F-A:]) 



(19) 



In the following we calculate the impact of an orientation depen- 
dent selection function on the galaxy bispectrum by introducing 
models for the anisotropic galaxy selection function which are 
based on symmetry considerations and motivated by the toy models 
of tidal alignment disc ussed in Sect. [3] First we extend the linear 
alignment model from iHiratal ( l2009h to second order in the den- 
sity field, and then construct a new model the anisotropic galaxy 
selection function due to quadratic alignment. 

3.1 Linear alignment 

In this subsection we construct a model for the anisotropic galaxy 
selection function 6 based on the assumptions that the large scale 
tidal fields induce a preferred direction in galaxy formation, and 
that the alignment is of linear order in the tidal field. Additionally 
we require the average of e(h\x) over the sky to vanish. Then the 
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only possible contraction of the tidal field with the viewing direc- 
tion h is 



e{h\x) 



AnGa^p^ia) 



Ain,«,(V,V,V -(5„)5(x), 



(20) 



where is the Newtonian potential, a is the scale factor, and where 
we have used the Poisson equation to write e in terms of the di- 
mensionless tidal field. A| is an expansion coefficient which en- 
codes the degree to which galaxy orientations are non-random due 
to tidal fields and the strength of galaxy orientation-dependent se- 
lection effects. Note that both effects need to be present in order to 
have Ai Q. 

To second order in the linear matter density field the 
anisotropic selection function in Fourier space can be written as 



mk) « A, 



(«.Af-ij[S"U) + 5<^)(A)]. 



(21) 



This expression is transformed to redshift space by Taylor expand- 
ing the real space expression and using Eqs. ([8l|9j 

e'Cnlx") = e{h\x) « 6(m|x') + (x - x') ■ VeiMx") + 0{6^) 

= e(h\x') + fh-W V-^6''\x'') h ■ Ve(«|jt:') , (22) 

and hence in Fourier space 

r<"(n|F) = e<'>(fj|F) 
r<^'(M|F) = e<2'(«|F) 



fd^ 



A:^)6<''(m|A:]). (23) 



Using this form for the selection function in combination with 
Eq. ( |19t , we now calculate the galaxy bispectrum modulated by 
linear tidal alignments. Then the first order observed density con- 
trast is given by 



5*<"(k^) = S<"(k) fc, - -A, + (A, + /);/! 



(24) 



The different terms contributing to the observed galaxy bispectrum 
can be calculated as 

»(kp5f <"(k^) (^)) = {27:f6^{ie^^^P{ie^)P{k'^ 



X (fci - -Ai +(Ai +f)n\ 



X \ ln - -A, +(Ai +f)n2 

X 2Z2(ki,k2) , 
"(k5)§f^<'>(k^) e»'''(nm)) = {2nf6^{k\2^)P(ie^)P{k\) 

1 



(25) 



X \b 



-Ai +(Ai +/)yu; 



X 1*1 - +(Ai + ,f)n2 



¥ I 1 

k\l 1 
+ AifiiiH2-ji:Ui2- 2 



convolution of first order density contrast and anisotropic selection 
function 

(5f <"(k^,)§f ^<"(k;) (Sf ® e-^'") (n, k^)) = {2nf6^ {k^J P{ie,)P{kl) 
X - ^Ai + (Ai + f )n]^ - ^Ai + (A, + f)^]^ 

X A, + [^12 - + (^1 + /y"^) - ^ j} . (27) 

Hence the galaxy bispectrum modulated by linear tidal alignments 
is given by 

Bf'\k\,kl k]) = [fo, - ^ + (Ai + /)//?] 
x[fei-^+(A, +/);/^j 

x|2Z2(Fi,A:^) + 2A, ^ - ^jF2(A:^A;^) 

+Ai foi J/i? + ^2 - ^) + I (6j"?/^2 - y"? - y"2) 

+^l,/'y"i^2 



^2/2 lU^W^ 1 



F(fci)P(yt2) + 2perm. 



(28) 



5. 7. 7 Transverse galaxy bispectrum 

As the full redshift space bispectrum is a complicated function of 
configurations described by 5 parameters (3 parameters specifying 
triangle shape, and 2 angles describing the orientation with respect 
to the line of sight), we will now simplify Eq. {28} by considering 
only triangles in the plane of the sky (/i, = 0), which are the easi- 
est to model and are the triangles observed in photometric redshift 
surveys. In this case, we find a galaxy bispectrum 



B\^Hk,,k2,k,) = {b, - ^) [2(fo, - -f)F2(ki,t 



+b2 --^^ib, 



P(ki)P{k2) + 2 peim. 



(29) 



Comparing this expression to Eq. ([6}, one finds that the effect of 
linear tidal alignments on the transverse galaxy bispectrum can be 
described as a rescaling of the galaxy bias parameters 



b\ , 

3 



(30) 



iHiratallioogI) found that the same rescaling of b\ applies to the real- 
space (/i, = 0) galaxy power spectrum. Therefore, the use of the 
real-space power spectrum and bispectrum to eliminate galaxy bias 
parameters and extract erg is robust against linear tidal alignments. 
However, this robustness does not extend to the pi + modes. 

For later use, we also write out the systematic error in the 
transverse galaxy bispectrum induced by linear alignment 



^B\^\kuk2,k,) 



2\b\A,-b, 



27 



FiikxM) 



-bs'^ 



P(ki)P(k2) +2 perm. 



(31) 



(26) 



3.1.2 Normalization 

FoUowing lffiratal j2009l) , we use Ai a -0.024 for LRG-type ellip- 
and the contribution from the last term in Eq. l |19b containing a tical galaxies. This is a rough estimate which is based on the as- 
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sumption that elliptical galaxies are optically thin triaxial systems, 
that the deviation from spherical symmetry can on average be re- 
lated to the tidal field (with correlation strength B), on different 
models for the orientation dependence of a galaxy's apparent mag- 
nitude (parametrized by;^'), and the slope of the galaxy luminosity 
function ?/: 



Ai=2tjxB. 



(32) 



While the total flux of an optically thin galaxy is not affected by 
tidal alignments, the average isophotal ellipticity and projected ef- 
fective radius of a galaxy become a function of the tidal field. 

The selection of galaxies in a survey will be modified by tidal 
alignment if part of the selection criteria is a magnitude cut, and if 
the apparent magnitude of a galaxy depends on its orientation. The 
apparent magnitude of a galaxy is nearly orientation independent if 
measured using Petrosian magnitudes or model magnitudes which 
are based on an accurate model for the radial profile, then at the 
level of the toy model considered by Hirata ( 2009) 6 » 0. 

If galaxies are selected using isophotal magnitudes or aperture 
magnitudes, more light will be counted if a galaxy is viewed along 
its long axis than its short axis. The selection factor in Eq. l |32t 
depends on the method used to measure galaxy fluxes (c.f. Fig. 2 
in|Hirata 2009), and it translates the fractional change in effective 
radius induced by intrinsic alignment to a fractional change in mea- 
sured flux. This change in measured flux moves galaxies across the 
selection threshold, and it is translated into change in number den- 
sity by assuming a luminosity function with slope -77. 

The strength of the tidal alignment efi'ect B is determined from 
measurements of t he density-ellipticity cross-correlation function 
jHirata et alj|2007h . Our chosen normalization further assumes a 
LRG luminosity function with 77 = 4.0 and galaxy selection based 
on isophotal magnitudes measured within ~ 3 effective radii. Also 
note that this normalization is based on observations around z = 0.3 
and should only be used near this redshift as the LRG luminosity 
function and the correlation between tidal field and galaxy orienta- 
tion may show strong evolution with redshift. 

3.2 Quadratic alignment 

The leading-order alignment of galactic angular momentum in tidal 
torque theories is quadratic in the tidal tensor because of the need 
for both a tidal field and an anisotropic inertia tensor on which it 
can act. 

The anisotropic selection function for a disc galaxy is gener- 
ally a function of its inclination ; (defined by cos i = L ■ h). While 
; is in the range < ; ^ ;r, we expect most selection criteria to 
be symmetric with respect to an observer being above or below the 
plane of the target, so it follows that the anisotropic part of the se- 
lection function contains only even-order spherical harmonics: 



r{Qh,x)= ^v/'j(cosO, 



(33) 



where is a Legendre polynomial. Using Eq. l ll7b . and noting 
that for a disk galaxy, we may replace the general integration over 
orientations Q €SO(3) with an integration over directions of the 
angular momentum vector L € 5^, we may write 



e(h\x) : 



2 - 1. 



p{L\x)Pjicoiii) A-L. 



(34) 



Because the quadratic alignment model contains two factors of the 
tidal field, which are spin 2, p(L\x) can contain spherical harmon- 
ics only through order 7 < 4. For simplicity, we will focus only on 



the quadrupolar J = 2 term in the sum (while noting that the hex- 
adecapolar alignment i = 4 is in principle possible). Then Eq. i34l 
implies that 



<E(h\x) oc {P2(L ■ h)). 



(35) 



where the average is taken over the local probability distribution of 
L. Equivalently, using Eq. dlSt . we find that 



(36) 



e{h\x) = A2 \hihj - -6ij^ Ti,,T,,j. 
We relate T, ; to the dimensionless shear field tensor T, ,, 



1 



1 ^kikj-^dijk'mk) 



kjkj 



(37) 



by approximating the scalar = TijTji with its expected value 
2 



= {T-) = -cr\R) 



(38) 



i.e. we approximate fjj a C^'Tij. As this expression for the 
anisotropic selection function is already second order in the den- 
sity field, effects associated with mapping 6 to redshift space only 
enter at higher orders than considered in this analysis and in the 
following we will drop the superscript s to denote Fourier modes in 
redshift space. 

Note that C~ is proportional to the variance of the smoothed 
density field smoothed on the halo collapse scale R, since the den- 
sity and tidal fields are both derived by taking second derivatives of 
the potential. 

Then the contribution of quadratic alignment to the orientation 
dependent selection function can be written as 



is. 



6^'\k')P\k"), 



(39) 



where k" = k - k' . This term contributes to the observed galaxy 
bispectrum via 

AB'^\k,,k2,k3) = 2A2 (*1 - y + (^1 +/)/'! ) 

xi^b,-^+(A,+f)^iijF(k,)P(k2) 

f , , 1 ^ » 2 

xl/Iifliki • A2 - - (/"i +y"2 + (^1 • *2)^) + q 

+ 2 perm. 



9 

(40) 



Here Aj t if the galaxy population under consideration is also 
subject to linear alignment, and we have defined A2 = A2/C^. 



3.2.1 Transverse galaxy bispectrum 

The quadratic alignment model modifies the observed transverse 
galaxy bispectrum by 



AB'^^-\kuk2,ki) = -A2b\ --[ki-k2f 

+ 2 perm. 



P(ki)P(k2) 



(41) 
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Figure 1. Effect of quadratic alignment on the reduced transverse galaxy 
bispectrum with h\ = \,ki = 0.05/! Mpc"', and where 612 denotes the 
angle between k\ and k2, and for A2 = 1. 
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Figure 2. Systematic offset of the reduced transverse galaxy bispectrum 
due to quadratic alignment with hi = \ and A2 = 1 as a function of triangle 
shape and scale. Shown are all possible closed triangle configurations with 
k\ ^ k2 > k3 for a given ki , areas in configuration space which do not 
correspond to a closed triangle are shown in white (located around the top 
and bottom left comer of each plot). Equilateral triangles ai'e located in the 
upper right corner of the configuration space, isosceles triangles lie on the 
upper diagonal, and collinear (9 [2 — » 0) triangles near the lower diagonal. 



Note that this systematic offset is independent of b2, and its ampli- 
tude scales linearly with At and quadratically with hi . The system- 
atic offset cannot be expressed as a simple rescaling of the galaxy 
bias parameters due to its shape dependence. Figure[T]illustrates its 
effect on the reduced transverse galaxy bispectrum 



Qg(ki,k2,k,)- 



BJki,k2,k3) 



PMi)PJk2) + PAki)P^(k3) + PAk2)PJh) ' 



(42) 

which is only mildly dependent on cosmology as the amplitude of 
fluctuations has been divided out. The shape and scale dependence 
of A2g is further illustrated in Fig.|2l which shows the systematic 
offset for all possible closed triangle configurations with ki'^ 
k3, with the left plot showing triangles with ki = 0.05/i/Mpc and 
the right plot showing triangles with ki = 0.2/i/Mpc. The system- 
atic offset is negative for triangles which are close to collinear, and 
for the scales considered in this analysis it shows little scale depen- 
dence. 



3.2.2 Normalization 

Similar to the normalization of the linear alignment contamination 
outlined in Section [3.1.2| the magnitude of the observed contamina- 
tion due to quadratic alignment again depends on (i) the orientation 
dependence of the recovered flux (continuum or line), (ii) the slope 
of the galaxy luminosity function, and (iii) the strength of the tidal 
alignment effect. We may use models for (i) and direct measure- 
ments for (ii), but (iii) is harder. For the linear alignment model 
we were able to use the observational constraints from the density- 
ellipticity cross-correlation function, but this is not an option here 
as the quadratic alignment contribution to two point statistics van- 
ish to leading order. Another option would be to set limits using 
the observed ellipticity variance, whi ch must set an upper li mit on 
(this was the approach followed in lCrittenden et al.l200ll for es- 
timating the intrinsic ellipticity coiTelation contamination of weak 
lensing surveys). We will take an even simpler approach here, and 
use some simple theoretical arguments on the value of a. 

In the tidal torque model, the distribution of disk normal 
vectors L given some tidal tensor T can be approximated by 
dCrittenden et alj|200ll) 



,j - 1 / 30- Qa- . - . 
p(L|T) « _ 1 + _ _ —L,LjnTjt 



(43) 



For a geometrically thin disk with normal vector L observed along 
the z axis, the inclination is cos ; = L3. The following constraints 
can be placed on a: 

• Since LiLjTjiiTj^ can take on any value between and |, the 
requirement that p(L|T) > sets the constraint |a| < |. 

• If one neglects correlations between the external tidal field 
and the mo ment of inertia ten sor of the collapsing protogalaxy, one 
finds tt = f jLee & Penll200(]|) . 

• The angular momentum of the disc of a galaxy may be dis- 
aligned from that of its host halo, due to e.g. torques between 
the disc and halo, or due to the disc containing only a specially 
selected subset of the halo's baryons. For a Gaussian distribu- 
tion of disalingment angles with rms per axis 0, the JM spher- 
ical harmonic component of p(L|T) is suppressed by a factor of 
exp[-y(y + 1)0-/2]; since we have a quadrupolar anisotropy {J = 
2), a is suppressed by a factor of exp(-30^). 

The above arguments suggest that \a\ of several tenths is plausible, 
but in no case should it exceed |. Also, while the simplest version 
of the tidal torque hypothesis implies a > 0, there is no physical 
reason why negative values should not be allowed. 

Next we determine the relation between an inclination depen- 
dent observed flux and the selection function e: Assume a galaxy 
flux distribution with slope dln/i/dlnFn,,,, = -rj. Then the number 
density of galaxies per logarithmic range in the intrinsic flux Fj per 
unit solid angle of disk orientation is 



M(F„L)ccFr"p(L\T). 



(44) 



Let the observed, inclination dependent flux be F(i) = FiO(i). The 
number density of galaxies above some threshold flux Fq then eval- 
uates to 

U{>Fo)cc fd'L r dlnFiFrV(ilt) 

J Jf„/(1>(,) 



Jo 



Fo/m 

mm" 



- 1 F2(COS(0) 



sin di, (45) 



where we have performed both the integral over and over Fj 
(since the latter is simply a power law), and defined Tj. = T^jTij. 
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Defining 



j^[<I>(0]'';'2(cos!)sinidi 



(46) 



£mi)f sin i d;- 
the anisotropic part of the observed galaxy count can be written as 

V 



9a 

T 



(47) 



Combining this with Eq. l l36t , we conclude that A2 = -r^aij/, and 
hence 



^ 9 a 27 aij/ 



(48) 



The top-hat variance is relate d to the bias of the galax ies if the mass 
function is nearly universal ( Sheth & Tormen|[l999l ); for example, 
at = 1 we have crjiR) = 2.96, whereas at = 2 we have crjiR) = 
0.83. 

The last step in obtaining a numerical estimate for A2 is eval- 
uating the orientation dependent selection factor ifi. This requires a 
model for the angular distribution of emitted radiance <t>{i), which 
also determines the selection probability p(i) oc [0(()]''. Several 
geometric toy models for the vertica l distributions of emitters and 
dust are discussed by iHiratal | |2009|) . and for galaxy distributions 
with 7} X 2 (appropriate for [On] and Hcf surveys), i/^ is found to 
be of order a few tenths: for example, it is i/r = 0.4 in the optically 
thick slab model; if/ = 0.23 (0.30) in the uniform slab model with 
normal optical depth r = 0.5 (1.0); and ij/ = 0.26 (0.37) in the 
sheet-in- slab model with r = 0.5 (1.0). 

These toy models suggest that A2 will be of order unity and we 
assume A2 = 1 for illustrative purpose in the following analysis^ 
For application to any survey the normalization must be calculated 
based on the detailed selection criteria and galaxy distribution. 



4 FISHER MATRIX ANALYSIS 

We now estimate the parameter bias induced by a tidal alignment 
contamination by performing a Fisher matrix analysis for a sur- 
vey with characteristics similar to the Dark Energy Survey (DESjl 
assuming that one would use the angular bispectrum of a slice of 
galaxies in photometric redshift space. A spectroscopic survey cov- 
ering a similar volume and oversampling the density field (nP > 1) 
would of course yield tighter constraints, but a full Fisher analysis 
of such a survey including redshift space distortions and finger-of- 
God parameters is beyond the scope of this paper 



4.1 Survey characteristics and analysis details 

Our fictitious survey has the same area as the DES, = 5000 
square degrees. We assume a constant comoving galaxy density 
over the redshift range of interest and use a radial galaxy sel ection 
function of the form expected for the DES jNock et ai]|2010h . 



dProb 

dz 



(49) 



In order to project out redshift space distortions we consider the an- 
gular clustering of galaxies projected over a finite radial distance. 

' In principle, either sign of A2 is allowed by our above calculations; for 
negative Ai the direction of the pai'ameter biases should be reversed. 
^ URL: http : //www . darkenergy survey . org/ 



For our theoretical modeling the projection over a finite range in 
radial distance is equivalent to a projection over a finite redshift 
range, and we choose 0.4 < z < 0.6. Observationally, this mapping 
is complicated by the distribution of photometric redshifts and the 
effect of redshift space disto rtions on the boundary of a region se- 
lected in redshift space (e.g. lPadmanabhan et al.ll2007l ; lNock et al.l 
|20l3). 



4.1.1 Binned angular multispectra and covariances 

We calculate the angular powe r and multispectra !Pjv using the Lim - 
ber equation in Fourier space ( lKaiseJl992l ; lBuchalter e"tai]|2000h : 



Hi ■ ■ ■ In) ■ 



r°\ r(x)„ Ih In . 
'^X -jjj^Pn\-,---,—,x] 

Jz=OA 



X 



X 



X 



(50) 



where is the three dimensional A'-point correlation function in 
Fourier space. In the following we use V, S, T to denote the 
angular galaxy power spectrum, bispectrum and trispectrum. 

For a linear alignment contamination, the change in the ob- 
served angular galaxy bispectrum is described by the same bias pa- 
rameter rescaling (Eq. J30t ) as for the transverse galaxy bispectrum 
discussed above. The magnitude of the systematic offset in the an- 
gular galaxy bispectrum induced by a quadratic alignment contami- 
nation is proportional to and independent of ^2- As the angular 
projection mixes different physical scales, the exact configuration 
dependence and normalization of the angular bispectrum contami- 
nat ion depends strongl y on the radial selection function (for details 
see lFrv& Thomas! 1999) . As can be seen from Fig.[2]the systematic 
olTset on the reduced transverse galaxy bispectrum is only weakly 
scale dependent, thus with our choice for the radial selection funt- 
cion the angular reduced bispectrum has very similar shape depen- 
dence. 

The Limber approximation requires the transverse scales un- 
der consideration to be significantly smaller than the radial pro- 
jection depth, hence we limit our analysis to angular scales cor- 
responding to comoving Fourier modes k ^ 0.04/; Mpc"'. As our 
intrinsic alignment toy models and biasing approximation are not 
designed to describe in the non-linear regime of structure forma- 
tion, we will only consider angular frequencies corresponding to 



0.04 hMpc-' < yt < 0.2 /iMpc" 



(51) 



We approximate the galaxy power spectrum by the linear matter 
power spectrum rescaled by the linear bias (Eq. |7); bispectra and 
trispectra on these scales are approximated by the tree-level pertur- 
bation theory in combination with local biasing (Eq. [5), i.e. using 
Eqs. Q, and l lA4t . These are evaluated using transfe r func- 
tions generated by CMBFAST dSeliak & Zaldarriagdll996[) for the 
best-fit WMAP 7 cosmology ("Komatsu et al."2010'). Compared to 
an approach combinin g the halo model with halo occupation distri- 
bution modeling (e.g. iBerlind & Weinberg I2OO2I ; ICoorav &Sheti3 
l2002i) this is computationally much faster, the only model input is 
our biasing prescription and does not require halo models for intrin- 
sic alignment. In the large scale limits the halo models asymptote to 
the perturbation theory result, and at the scales of our analysis the 
galax y power spectrum is fairly well described by pert urbation the- 
ory ( Cooravl2004ISmith et al.l2008l) . At redshift z = O. ISmith et all 
(200^ find the reduced halo model bispectrum with k2 - 2ki to 
be in close agreement with perturbation theory results at scales 
ki < 0. 1 /z/Mpc, except for collinear configuration {On — » 0). As 
we only consider triangle configurations with all angular frequen- 
cies k\xi < 0.2 /jMpc"', the perturbation theory results should be 
sufficient at the level of this analysis. However, at scales smaller 
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than/t ~ Q.ih /Mpc ISmith et al J j20o3) and lOuo & Ji"ri3 i2009l) find 
the bispectrum measured from simulations to differ at the 10-20% 
level from the perturbation theory. Note that these systematic ef- 
fects on the determination of bias parameters on small scales are 
larger than the tidal alignment contaminations discussed here. 

We model the observed power spectrum by averaging the an- 
gular power spectrum over bins of width A/, 

W+1/2A/ I 
Jl-i/lAI lCs.1 

and the corresponding covariance is given by 

n2 



(52) 



Q. 











n 



(53) 
(54) 



where h is the average projected density of the galaxy population 
under consideration. Here the first term is a combination of Gaus- 
sian cosmic variance and shot noise. The second term involving 
the trispectrum of parallelogram configurations is the non-Gaussian 
power spectrum covariance. 

The bispectrum is sampled with uniform binning AZ in all an- 
gular frequencies. Defining 

U (55) 

Ji J/,-l/2A/ Z;A/ 

the bin-averaged bispectrum is given by 

Sduhh) = ^JJmuh,h)6o{h +h + hi (56) 

We approximate the expression from ljoachimi et aL I( l2009l) for the 
full non-Gaussian covariance of the bin-averaged bispectrum by 



X Di. 



1 



nh) + - 

n 



Pill) + T 
n 



nh) + r 

n 



■'li.h.h.hJs.k 

S2 ' Jl J2 JT, J5 J6 

X Iso (h + h + i6)mj2, hmh, h,k) 



ShM [ f f f [Soih 

Ji J2 Ji Js Je 



+ Soi-h+h+hJTduhJs, 16)^13)} + ^Veim. 



+ I2 + h) 



(57) 



where the symmetry factor £);- is non-zero only for diagonal el- 
ements of the covariance {[hJiJi) - {UJsJe])'- Aj.../,, = 1' 2, 
or 6 for scalene, isosceles, or equilateral triangles respectively. If 
h,l2, h form a triangle, then A"'(Zi, Z2, h) is the area of this trian- 
gle, otherwise A"' = 0. The first term is the Gaussian (diagonal) 
part of the covariance which is proportional to the product of three 
power spectra which have been modified to account for Gaussian 
shot noise. The second/ third terms are non-Gaussian contributions 
from triangle pairs which have at least one common side so that the 
pentaspectrum can be factorized into two bispectra/ a trispectrum 
and a power spectrum. We have dropped a term which is propor- 
tional to the general connected pentaspectrum. 

4.2 Biased parameter estimates for galaxy bias parameters 

Having set up a model for the observable data and their covari- 
ances, we can now quantify the power of our fictitious survey at 



constraining model parameters using the Fisher matrix 
r„, = |^Cov-'(A^)^.|^Cov-'(^:s)|^, 



(58) 



where the P and S are data vectors with the binned angular galaxy 
power spectrum and bispectrum as data points. The data vectors and 
their covariances depend explicitly on the bias parameters through 
Eqs. l?] (6] IA4t . Note that we do not include cross-correlations 
between power spectrum and bispectrum, both for simplicity and 
because they are sm all in the weakly nonlinear regime (but see 
ISefusatti et alj2006l for their constraining power in the weakly non- 
linear regime). The parameters of interest here are the linear and 
quadratic galaxy bias and we marginalize over the normalization 
of the matter power spectrum erg, i.e. p = {bi,b2,(T%). Our fiducial 
model assumes erg = 0.8, no intrinsic alignment contamination, 
and covers a range of bias parameters, while all other cosmological 
parameters are fixed to their best-fit WMAP 7 values. 

The inverse Fisher matrix serves as a lower limit on the 
marginalized covariance of statistical parameter errors 



{6p„6pp) = [t-') 



(59) 



Hence the statistical error on the inferred parameters is inversely 
proportional to Vf2, as can be seen from the expressions (Eqs. 1541 
|57t for the data covariances. The presence of a systematic error AS, 
AV in the data which is not included in the model induces a bias in 
the parameter estimate compared to its fiducial values. To first order 
it is given by (e.g. ,Huterer et al. 2006; Arnara & Refregier 2008.) 

Apa = {p„} - pI' = aP' Gov-' {'P, r) — 



A'B' Cov^' 



dS 
dpp 



(60) 



where the data vectors and covariances are evaluated at the fiducial 
model. 

This systematic bias is independent of the survey area, but it 
is influenced by our choice of survey parameters through the se- 
lection function (Eq.|49t and data binning scheme. It also depends 
on projected number density of the galaxy population of interest 
as h determines the importance of shot noise. We adopt a uni- 
form sampling with 20 equidistant bins in all angular frequencies 
(h,l2,h) corresponding to Eq. l lSlt and assume a projected den- 
sity of n = 1 /arcmin^ for a galaxy population in the redshift range 
0.4 ^ z < 0.6. 

The systematic error on the bispectrum, AS, due to linear 
or quadratic alignment is modeled by the line of sight projection 
(Ea.l50t of the tidal alignment contaminations (Eqs. |29ll4U calcu- 
lated in Sect.|3] We set AP = for the quadratic alignment model 
as the first correction to the power spectrum is third order in the 
density contrast. In agreement with our findings from Eq. OOt . the 
systematic error induced by linear ali gnment on the galaxy power 
spectrum is given by (cf. lHiratall2009h 



(61) 



where we have restricted k to be orthogonal to the line of sight as 
only these modes survive the Limber approximation. 

Figure [3] shows the marginalized Fisher matrix estimates of 
statistical parameter errors (95% C.L.) obtained with our fictitious 
survey in the absence of an intrinsic alignment, and the systematic 
bias induced by a linear or quadratic alignment contamination. 
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Figure 3. Systematic errors induced by intrinsic alignment. Ellipses show 
95% C.L. statistical errors on parameter estimates in a DES-like surveys 
for a fiducial model with erg = 0.8, for a galaxy population with h[ = 
1 (top panels) or =2 (bottom panel) and bi e {-0.5,0,0.5). Open/ 
filled arrows illustrate the systematic parameter shift induced by a quadratic/ 
linear intrinsic aUgimient contamination. 
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Figure 4. Removal of quadratic alignment bias. Dotted ellipses show the 
biased parameter estimates and their 95% contour regions in the presence 
of quadratic alignment contamination with A2 = 1 which is unaccounted 
for in the analysis. The solid ellipses illustrate the 95% contour regions of 
the unbiased parameter estimates in an analysis which includes a quadratic 
alignment contamination and marginalizes over A2 . 



The systematic bias induced by a linear alignment contami- 
nation (solid arrows) we find through the Fisher matrix analysis 
(Eq.|60ll is in agreement with the analytic result (Eq.|30t. The pa- 
rameter bias on bi is independent of the value of ^2 assumed in the 
fiducial model and the solid arrows of different color are indistin- 
guishable. Assuming a normalization of A] = -0.024 as discussed 
in Sect. [3] the systematic error on ^2 is comparable to the 95% 
C.L. statistical error for b2 in our survey. The systematic error on 
hi caused by the linear alignment model is smaller, but may still be 
important if many photo-z slices are used in the parameter analysis. 
In the limit of our toy model, the effect of linear alignment on the 
angular galaxy power spectrum and bispectrum is fully described 
by a systematic error in the linear and non-linear bias parameter 
(Eq.|30t and it has no effect on measurements of ag. 

The strength of the quadratic alignment contamination de- 
pends on triangle shape and size; it is not well described by a 
rescaling of the galaxy bias parameters. Hence the Fisher matrix 
estimates for the systematic parameter errors depend on the bin- 
ning scheme and range of scales adopted in the analysis. For our 
choice of 20 equidistant bins per angular frequency, and with the 
range of scales of 0.04-0.2/! Mpc"', we a systematic shift towards 
larger non-linear bias ^2 and smaller bi. The latter is degenerate 
between bi and Cg. The plot illustrates a quadratic alignment con- 
tamination with normalization Ai = 1. As can been seen from 
Eqs. I I41I |60t , the systematic bias is linear in A2, and it reverses 
sign if A2 < 0. While exact form of the systematic error caused by 
the toy model for quadratic alignment depends on a number of pa- 
rameters, it may cause a significant contamination in our fictitious 
survey if IA2I > 0.5, or if (as we expect) multiple photo-z slices are 
used to reduce statistical errors. 



5 DISCUSSION 

Using simple toy models for intrinsic alignment and the local bias 
approximation we have analyzed the effect of tidal alignment on 
the galaxy bispectrum. If the orientation of galaxies depends on the 
surrounding tidal field, and if the detection probability for galax- 
ies is orientation dependent, the observed clustering of galaxies is 
modified by tidal alignments. This astrophysical contaminant can 
introduce systematic errors to parameters derived from the bispec- 
trum. 

A toy model for linear aligiraients jCatelan et al.l200ll) . which 
is based on the assumption that tidal fields elongate/compress 
haloes and thus determine galaxy shapes, results in a rescaling of 
linear and non-linear galaxy bias parameters that is proportional to 
the strength of the halo shape distortion. The presence of this sys- 
tematic error in the observed galaxy bias measurements cannot be 
detected from projected clustering data as the strength of the align- 
ment contamination is completely degenerate with the unobserv- 
able true bias parameters and outside information will be necessary 
to remove it. Normalizing the strength of the linear tidal alignment 
toy model to measurements of intrinsic alignments in weak lensing 
observations, we find that linear alignment may introduce system- 
atic errors to galaxy bias measurements at the percent level (again 
using only the real-space observables), and thus will likely not be 
significant. 

Using a simple model for quadratic alignment based on galaxy 
spin correlations in linear tidal torque theory we calculate a system- 
atic contamination which modifies the shape of the galaxy bispec- 
trum. Depending on survey characteristics, we find that quadratic 
alignment may introduce significant systematic errors to the galaxy 
bias parameters and the normalization of the power spectrum de- 
rived from the angular galaxy bispectrum. As the quadratic align- 
ment contamination has different shape than the galaxy bispectrum, 
one can include a model for the contamination in the analysis and 
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marginalize over its normalization. Figure |4] illustrates how such a 
marginalization may remove the systematic bias at the cost of larger 
statistical errors. The biased data points and contour levels (dashed 
lines) are taken from Fig. [3] for a fiducial model with b[ = I and 
b2 = 0. The new statistical errors including marginalization over A2 
are calculated by adding A2 as a nuisance parameter and including 
the contamination signal in the fiducial model of the Fisher matrix 
analysis (S S + AS in Eq.l58]l. 

This analysis lives in the weakly nonlinear regime to enable 
the use of simple models for linear and quadratic alignment. As the 
information content of the bispectrum increases dramatically with 
the maximal spatial frequency that is included in an analysis, any 
realistic analysis will have extend well into the quasilinear regime. 
While mode ls from the redshift space bispectrum on these scales 
jSmith et alj|200 8) approach the required accuracy for such anal- 
yses, the treatment of tidal alignments including the non-Gaussian 
nature of the angular moment distribution and non-linear stages of 
galaxy formation requires further work. 
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APPENDIX A: TREE-LEVEL GALAXY TRISPECTRUM 

To calculate the tree-level matter trispectrum we need to con- 
sider the density contrast to third order as the tree-level Trispec- 
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trum splits into two types of connected terms, ^S*''(5*'*5'^'5*^'^ and 
/^(i)g(i)^(i)g(3)^ The third order density contrast is given by jprvl 
Ii984h 

xWKk^W\k2)W\k -ki- ki) , (Al) 

with the third order coupling function Fj . One finds for the matter 
trispectrum 

(27r)'5D(*1234)7'pt(*l,*2,*3, A^4) ~ 

{6^'\k^)d^'\k2W\ki)S^^\k^)) + 3 perm. 
+ {s'^'\ki)6^'\k2)S'^^\ki)P\k4)) + 5 perm. (A2) 
After some algebra one obtains 

Tpt{ki,k2, ki, k^) = 6Fl(ki,k2, ki)P{h)P(k2)P{h) + 3 perm. 

+A[P{kn)F2{ku-kn)F2(k2,ka) 

+P(k2i)F2(kuk23)F2(k2-k23)] 

xP(k,)P(k2) + 5 penn. (A3) 

If one assume the third order galaxy bias (^3) to be zero, two 
types of additional terms containing the quadratic galaxy bias con- 
tribute to the galaxy trispectrum, fc^S*" 

and (foiS'" foiS'" fc2§"' ® S<" fc2S<" Hence our model for 

the galaxy trispectrum is given by 

rgal(*l,A2,*3,*4) ~ b^iTj„(k,,k2,ki,k4) 

+2b]b2P{h)P{k2)[P{hi)F2iku-kn) 
+P(k24)F2{k2,-k23)] + 5 perm. 

+Ab\blP(h)P{k2) [P{kn) + P{k23)] 

+5 perm. (A4) 
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